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FOREWORD 


This  study  was  conducted  for  the  Naval  Air  Development 
Center,  Warminster,  Pennsylvania  by  the  Aerodynamics  Department  of 
the  McDonnell  Aircraft  Company  (MCAIR) ,  St.  Louis,  Missouri.  The 
objectives  of  this  study  were  to  determine  the  fountain  impinge¬ 
ment  momentum  flux  transfer  coefficients  associated  with  three 
V/STOL  aircraft  hovering  at  several  aircraft  heights  in  ground 
effect.  The  three  V/STOL  aircraft  considered  were:  (1)  the  MCAIR 
YAV-8B,  (2)  the  MCAIR  Model  260  and  (3)  the  Grunman  Model  698. 

The  fountain  impingement  momentum  flux  transfer  coefficients 
were  calculated  indirectly  by  applying  the  MCAIR  v/STOL  force  and 
moment  prediction  methodology  which  separates  the  suckdown  effects 
from  the  fountain  impingement  effects.  This  methodology  calcu¬ 
lates  the  suckdown  force  utilizing  a  potential  flow  panel  method. 
The  fountain  impingement  force  is  calculated  using  a  semi- 
empirical  methodology  that  calculates  free  jet  and  wall  jet 
properties  (including  jet  entrainment  velocities),  predicts  the 
ground  wall  jet  stagnation  line  location  between  impinging  jet 
pairs  and  calculates  2-jet  fountain  upwash  inclination  and  upwash 
momentum  flux  magnitudes.  The  calculated  suckdown  and  fountain 
impingement  forces  were  then  compared  with  the  total  induced  lift 
determined  from  wind  tunnel  tests  to  "back  out"  the  fountain 
impingement  momentum  flux  transfer  coefficient  for  each  of  the 
three  aircraft. 

The  principal  investigator  of  this  study  was 
Mr.  Lloyd  W.  Glaze  of  the  McDonnell  Aircraft  Company  Aerodynamics 
Department.  The  Naval  Air  Development  Center  contract  monitor  was 
Mr.  Marvin  M.  Walters. 

The  authors  gratefully  acknowledge  the  NASA  Ames  Research 
Center  for  their  cooperation  in  the  conduct  of  this  study  and  the 
Grumman  Aircraft  Corporation  for  providing  configuration 
information  and  test  results  for  the  Model  698  V/STOL  aircraft. 
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SECTION  I 
INTRODUCTION 


The  successful  design  and  development  of  V/STOL  aircraft 
equipped  with  powered  jet  lift  systems  requires  a  knowledge  of 
the  complex  flow  field  interactions  produced  by  these  multiple  jet 
lift  systems  in  the  hover  mode  both  in  and  out  of  ground  effect 
(Figure  1).  The  accurate  prediction  of  the  induced  forces  and 
moments  on  the  aircraft  is  necessary  in  order  to  assure  that  the 
propulsion  system,  control  system,  and  other  subsystems  may  be 
properly  designed.  A  key  element  in  the  analytical  modeling  of 
these  flow  fields  in  ground  effect  is  the  formation  and  develop¬ 
ment  of  fountain  upwash  flows .  The  prediction  of  the  trajectory 
of  the  fountain  upwash  is  also  of  importance  for  the  computation 
of  moments  and  in  order  that  lift  jet  system  exhaust  reingestion 
may  be  minimized.  A  considerable  amount  of  progress  has  been  made 
by  the  McDonnell  Aircraft  Company  (MCAIR)  toward  the  prediction  of 
the  complete  forces  and  moments  induced  on  a  V/STOL  aircraft  oper¬ 
ating  in  the  hover  mode  through  dedicated  experimental  and 
analytical  studies  conducted  through  various  Contracted  Research 
and  Development  (CRAD),  as  well  as  MCAIR  Independent  Research  and 
Development  (IRAD)  programs. 


amom-M 

Figure  1.  Flowfleld  About  a  VTOL  Aircraft  Hovering  in  Ground  Effect 


These  research  efforts  resulted  in  the  development  and  refine¬ 
ment  of  the  MCAIR  multi-jet  V/STOL  aircraft  force  and  moment 
methodology,  developed  for  the  Naval  Air  Development  Center 
(NAVAIRDEVCEN)  in  Reference  (1).  Originally  developed  for  round 
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jets,  this  methodology  was  subsequently  expanded  under  contracts 
to  the  NASA  Ames  Research  Center  to  include  propulsive  lift  sys¬ 
tems  employing  rectangular  nozzles  with  exit  area  aspect  ratios 
between  one  and  eight  (References  2  and  3).  The  approach  in  this 
methodology  is  to  separate  the  jet  entrainment  induced  suckdown 
forces  from  the  fountain  impingement  forces.  The  former  are  deter¬ 
mined  from  the  jet  entrainment  induced  potential  flow  field  util¬ 
izing  a  multiple  boundary  condition  panel  methodology,  while  the 
latter  are  determined  from  empirical  jet  impingement  data  and 
momentum  conservation  principles  with  the  inclusion  of  empirically 
determined  fountain  upwash  formation  momentum  loss  and  fountain 
upwash  airframe  impingement  coefficients.  This  empirical/theoreti- 
cal  analysis  is  embodied  in  the  MCAIR  Aircraft/Ground/Stagnation 
Line/Fountain  (AGSF)  methodology  (Reference  4). 

A  recent  experimental  investigation  conducted  by  MCAIR  for 
the  Office  of  Naval  Research  (ONR)  (Reference  5)  yielded  two 
fundamentally  significant  modeling  elements:  (1)  experimentally 
determined  fountain  formation  momentum  flux  recovery  coefficients, 
and  (2)  a  verified  analytical  model  for  the  fountain  upwash  flow 
trajectory.  The  fountain  momentum  flux  recovery  coefficient,  Aft' 
is  the  ratio  of  local  momentum  flux  in  the  fountain  upwash  to  the 
total  local  momentum  flux  entering  the  fountain  base  from  the 
ground  wall  jets.  A#  is  shown  as  a  function  of  jet  pair  exit 
momentum  flux  ratio  (or  thrust  bias),  MjeLQy/MjejnGjj  in  Figure  2. 

These  data  indicate  that  approximately  45%  of  the  momentum  flux 
entering  the  fountain  formation  region  between  two  equal,  verti¬ 
cally  impinging  jets  is  not  available  in  the  fountain  upwash  flow, 
(A^  s  0.55).  Figure  3  presents  the  variation  of  the  2-jet  foun¬ 
tain  inclination  (w  )  as  a  function  of  the  jet  pair  exit  momentum 
flux  ratio  for  various  values  of  A^  along  with  experimental  verifi¬ 
cation.  The  dashed  line  in  Figure  3  represents  the  fountain 
inclination  based  on  the  curve  fit  value  of  A^  shown  in  Figure  2. 

The  establishment  of  the  above  ground  flow  field  modeling 
elements,  together  with  recent  improvements  made  to  potential  flow 
panel  methods  for  the  determination  of  jet  entrainment  induced 
suckdown,  leaves  only  one  empirical  element  remaining  to  coiqplete 
the  prediction  methodology  for  a  wide  range  of  V/STOL  aircraft 
configurations.  This  remaining  element  is  the  quantification  of 
the  amount  of  momentum  flux  transferred  to  the  airframe  under¬ 
surface  from  the  upwash  momentum  flux  in  the  fountain.  This 
quantity,  Afj,  is  termed  the  fountain  upwash  momentum  flux  trans¬ 
fer  coefficient. 

The  fountain  impingement  momentum  flux  transfer  coefficient 
was  determined  for  the  YAV-8B  in  hover  at  one  height  in  ground 
effect  as  part  of  a  recent  MCAIR  Independent  Research  and  Develop¬ 
ment  (IRAD)  Project.  The  present  contract  study  includes  the 
calculation  of  A  f j  for  the  YAV-8B  aircraft  at  three  additional 
heights  in  ground  effect.  Furthermore,  Afj  was  calculated  for  the 
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SECTION  II 
APPROACH 


The  MCAIR  V/STOL  hover  methodology  is  illustrated  for  a  two- 
jet  configuration  in  Figure  5.  The  total  thrust  non-dimens ional- 
ized  jet  induced  lift  (Al/t)  is  defined  as  the  sum  of  the  non¬ 
dimens  ionali  zed  jet  entrainment  induced  suckdown  force  (Fs/T)  and 
the  non-dimensionalized  fountain  impingement  force  (Ffi/T)«  The 

fountain  impingement  force  is  in  turn  defined  as  the  product  of 
the  fountain  impingement  momentum  flux  transfer  coefficient  (X^j) 
and  the  total  fountain  upwash  momentum  flux  (Afi/T)  incident  upon 
the  aircraft;  thus. 


AL 

T 


Fs  ,  flfl 

ip  ^  f  I 


(1) 


The  approach  taken  in  this  study  is  to  calculate  the  suckdown 
force,  Fg/T,  and  the  fountain  upwash  momentum  flux,  and  by 

applying  equation  (1),  compare  their  sum  to  experimentally  deter¬ 
mined  values  of  the  total  induced  lift,  AL/t,  to  “back  out“  the 
corresponding  value  of  X^j.  This  was  done  for  the  three  aircraft 
illustrated  in  Figure  4  for  at  least  three  aircraft  heights  above 
ground  each.  The  procedures  used  to  calculate  Fg/  and  ftfj/  are 
described  below.  T  T 


1.  SUCKDOWN  FORCE  CALCULATIONS  -  To  calculate  the  suckdown  force 
on  a  V/S*fOL  aircraft  hovering  in  ground  effect,  the  aircraft,  lift 
jet  and  wall  jet  surfaces  are  paneled  for  use  with  a  surface  singu¬ 
larity  panel  method.  Figure  6  presents  the  paneled  representation 
of  the  YAV-8B  in  hover  in  ground  effect.  The  panel  methodology 
selected  for  use  in  this  study  is  a  modified  version  of  the  PAN 
AIR  surface  singularity  methodology  of  Reference  6.  The  modifica¬ 
tion  permits  the  specification  of  non-zero  normal  velocity  bound¬ 
ary  conditions  on  jet  surface  panels  to  represent  jet  entrainment 
flow.  Zero  flow  normal  to  the  panel  surface  is  prescribed  as  the 
boundary  condition  on  aircraft  surfaces .  The  suckdown  methodology 
was  verified  and  the  entrainment  models  utilized  by  the  AGSF  pro¬ 
gram  were  investigated  prior  to  applying  the  methodology  to  the 
complete  aircraft  configurations. 

a.  Verification  of  Suckdown  Methodology  -  In  order  to  verify 
the  validity  of  this  approach  when  applied  to  aircraft  in  hover 
(V^O),  a  simple  test  case  consisting  of  an  axisymmetric  jet 
issuing  from  an  infinite  flat  plate  into  static  air  was  modeled. 
The  pressure  distribution  on  the  surface  of  the  plate  was  calcul¬ 
ated  using  two  different,  empirically  determined,  axisymmetric  jet 
mass  entrainment  models.  Figure  7  presents  the  entrained  mass 
flow  distributions  determined  for  a  circular  jet  in  five  separate 
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previous  investigations.  The  wide  disparity  in  the  curves  illus¬ 
trates  the  uncertainty  in  the  measurement  of  such  data.  Figure  8 
compares  the  pressure  distributions  calculated  by  a  panel  method 
using  the  entrainment  curves  of  wygnanski  and  of  Kleis  and  Foss 
with  experimental  plate  surface  pressure  data  determined  by 
Wygnanski.  The  curves  illustrate  the  effect  of  varying  the  free 
jet  mass  entrainment  rate  on  the  flat  plate  pressure  distribution. 


AIRCRAPT/QROUNO/STAQNATION  UNEJfOUNTAIN 
(AQSF)  PROGRAM  CALCULATES S 

•  FREE  JET  ENTRAPMENT  VELOCITIES 
e  WALL  JET  ENTRAINMENT  VELOCITIES 


#  JET  IMPINGEMENT  POINTS 

e  STAGNATION  LINE  LOCATIONS 
e  FOUNTAIN  UPWASH  DIRECTION 

•  FOUNTAIN  UPWASH  MOMENTUM  FLUX 
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Figure  5.  MCAIR  V/STOL  Ho vor  Methodology 
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Figure  6.  Panelsd  Representation  of  YAV-8B  in  Hover  In  Ground  Effect 
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Figure  7.  Entrained  Maas  Flow  Distribution  for  Circular  Jats 
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In  both  of  the  calculated  cases  shown  in  Figure  8,  the  mass 
entrainment  rate  was  considered  to  be  constant  within  the  poten¬ 
tial  core  region  (Z/D  ^  6)  and  constant  outside  the  potential  core 
region  (Z/D  >  6)  with  the  higher  entrainment  rate  occurring  in  the 
fully  developed  region.  The  accuracy  with  which  the  panel  method 
(in  conjunction  with  Wygnanski's  entrainment  model)  reproduces 
Wygnanski's  empirical  curve  demonstrates  the  validity  of  the  panel 
method  for  use  in  hover. 


Figure  8.  Pressure  Distribution  Due  to  an  Axlsymmetrte  Jet  Issuing  from  s  Wsll 

The  calculated  pressure  distributions  were  integrated  as  a 
function  of  radial  distance  along  the  plate  to  yield  a  normal 
force  coefficient  (Cjj) .  The  variation  of  integrated  normal  force 
with  radial  distance  is  illustrated  in  Figure  9  for  the  two 
entrainment  models.  Also  shown  in  Figure  9  is  the  normal  force  on 
the  plate  induced  by  a  single  circular  jet  (operating  at  the  same 
conditions  as  the  test  case)  determined  from  the  empirical  rela¬ 
tion  derived  by  Byrne  (Reference  12).  This  relationship  (shown  in 
Figure  10)  represents  a  correlation  of  induced  lift  data  obtained 
experimentally  for  several  V/STOL  aircraft  models  in  hover  out  of 
ground  effect.  The  normal  force  curves  of  Figure  9  indicate  that 
both  entrainment  models  tend  to  underpredict  the  actual  jet 
induced  force.  It  was  concluded  from  this  study  that  the  poten¬ 
tial  flow  representation  of  jet  induced  suckdown  is  valid,  but  the 
existing  empirical  entrainment  distributions  are  too  weak  and  must 
be  corrected  to  accurately  model  jet  induced  forces.  As  a  result 
of  this  finding,  an  experimental  test  program  was  conducted  using 
a  new  technique  to  indirectly  measure  axisynmetric  free  jet 
entrainment  using  a  hot  film  amemcmeter.  The  experimental  test 
program  and  test  results  are  discussed  in  Appendix  A.  A  procedure 
for  correcting  existing  entrainment  distributions  for  calculating 
suckdown  forces  in  hover  in  ground  effect  (IGE)  is  presented  in 
Section  II.l.C. 
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Figure  9.  Normal  Force  Oua  to  an  Axiaymmatric 
Jot  laauing  from  a  Flat  Plata 


Figure  10.  Hover  Lift  Loaa  Correlation 
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b.  Effect  of  Lift  Jet  Paneling  -  The  effect  of  the  lift  jet 
paneling  geometry  was  investigated  using  the  paneled  model  of  the 
YAV-8B  aircraft  in  hover  out  of  ground  effect  (OGE) .  Figure  11 
illustrates  the  paneled  aircraft  with  two  possible  paneled  jet 
models.  In  one  case,  the  paneling  represents  the  spreading  jet 
geometry  calculated  by  the  AGSF  program.  In  the  other  case,  an 
equivalent  cylindrical  jet  paneling  was  applied.  The  entrainment 
data  of  Kleis  and  Foss  was  utilized  for  both  cases.  It  was  veri¬ 
fied  that  both  paneled  jet  representations  generate  the  same 

calculated  suckdown  force  on  the  aircraft,  •”(  *  -.006.  This 

confirms  that  the  simple  cylindrical  jet  paneling  is  entirely 
adequate,  provided  that  the  jet  normal  velocity  (entrainment) 
distribution  is  adjusted  to  account  for  the  smaller  jet  perimeter. 
However,  the  experimentally  measured  suckdown  force  on  the  YAV-8B 
in  hover  OGE  is  -.020,  substantially  larger  than  predicted.  This 
result  reemphasizes  the  result  determined  using  the  infinite  flat 
plate  model  that  existing  empirical  free  jet  entrainment  rates 
underestimate  actual  entrainment  and  require  adjustment. 


DETAILED  JET  MODEL  SIMPLE  JET  MODEL 


Figure  11.  YAV-8B  in  Hover  Out  of  Ground  Effect 

Alternate  Jet  Models 


NADC-81106-60 


c .  Corrected  Free  Jet  Entrainment  Velocities  -  Due  to  the 
uncertainty  In  empirical  distributions  of  free  Jet  entrainment 
rates  (Figure  7),  the  free  jet  entrainment  velocities  utilized  to 
calculate  suckdown  forces  in  hover  IGE  are  based  on  corrected 
values  of  the  free  jet  entrainment  velocities  (Vfen^)  utilized  for 
hover  OGE.  In  this  procedure,  the  suckdown  force  out  of  ground 
effect  (  a)  is  calculated  using  the  entrainment  data  of  Kleis  and 

Foss.  Tne  result  is  then  compared  with  experimental  data  (or,  if 
such  data  is  not  readily  available,  with  the  general  curve  devel¬ 
oped  by  Byrne  [Reference  12],  Figure  10)  to  yield  the  corrected 
free  jet  entrainment  velocity,  V'fent,  as  follows: 

V'fent  Vf6nt  (°Ut  °f  9rOUnd  effect)  (2) 


d.  Evaluation  of  Wall  Jet  Entrainment  Velocities  -  As  origin¬ 
ally  calculated  'Ey  the  AGSF  program,  the  wall  jet  entrainment 
velocities  (vwept^  varied  directly  with  the  radius  of  the  jet 
impingement  region  (Rj)  and  the  free  jet  centerline  velocity  at 
impingement  (Vjcl)  and  inversely  with  radius  (r);  i.e.. 


vwent 


a  ( 


RIVjCl) 


(3) 


Although  this  relation  provides  the  correct  variation  of  Vwent 
with  radial  distance  (Figure  12),  the  variation  of  Vwerrt  with 
nozzle  height  (h/D)  is  the  opposite  of  that  which  would  be 
expected.  This  unreasonable  variation  was  found  to  result  from 
the  fact  that  Rj  increases  more  rapidly  than  Vjcj  decreases  with 
respect  to  increasing  h/D.  The  expression  for  VjcI  is  based  upon 
empirical  data  and  is  generally  assumed  to  be  accurate.  Conse¬ 
quently,  Vw^nt  was  redefined,  in  a  manner  consistent  with  the 
formulation  in  Reference  1,  as: 


vwent 


a  (. 


2.4D 


Ii£l) 


(4) 


(where  2.4D  represents  the  radius  of  the  impingement  region  at 
h/D  «  4).  Figure  13  illustrates  that  the  revised  wall  jet 
entrainment  velocities  exhibit  the  correct  variation  of  Vwen^  with 
h/D.  It  should  also  be  noted  that  Vwent  is  invariant  with  h/D  for 
h/D  <_  6.  The  reason  for  this  is  that  the  impinging  jet  centerline 
flow  is  still  within  the  potential  core,  where  Vjcj  is  constant. 
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In  order  to  verify  the  accuracy  of  the  magnitude  of  Vwent<  a 
single  test  case  was  devised  for  use  with  the  panel  methodology. 
The  test  case  (shown  schematically  in  Figure  14)  consisted  of  a 
circular  disk  (with  diameter  6)  emitting  an  axisymmetric  jet  (with 
diameter  D)  in  ground  effect.  The  free  jet  entrainment  velocities 
used  were  the  corrected  free  jet  entrainment  velocities,  V'fent, 
defined  by  equation  (2).  The  suckdown  force  on  the  disk  was  then 
calculated  using  the  wall  jet  entrainment  velocities,  Vwent,  of 
Figure  13  multiplied  by  various  constant  scaling  factors, 

The  variation  of  the  calculated  induced  lift  Fs/T,  with  is 
shown  in  Figure  15  for  a  disk-to-jet  diameter  ratio  (5/D)  of  8  and 
for  disk  heights  of  two  and  four  jet  diameters  above  the  ground 
plane.  The  dashed  lines  in  the  figure  represent  the  induced  lift 
determined  using  the  following  empirical  relation  (from  Reference 
12): 


J*  g  f  j  .M 

~  =  -  .015  [ - =- - £__]  exp  {-  [2.2  -  .24  (NPR-1)]  }  (5) 

(D/De-1) 

The  diameter,  De,  represents  the  equivalent  single  nozzle  exit 
diameter  and  is  defined  as: 


D 


(6) 


where  Ajei  is  the  nozzle  exit  area  of  the  i  th  nozzle  and  n  is  the 
number  of  jets.  It  should  be  noted  that  for  the  single  axisym¬ 
metric  nozzle  shown  in  Figure  14,  D  is  equivalent  to  De. 


The  points  in  Figure  15  at  which  the  curves  cross  the  dashed 
lines  determine  the  value  of  required  for  adjustment  of  the 
uncorrected  wall  jet  entrainment  velocities  (Vwent)  to  yield  the 
empirical  value  of  the  induced  lift.  Curves  such  as  those  shown 
in  Figure  15  were  also  obtained  for  D/De  =  6  and  D/De  =  12  at 
h/De  =  2  and  4.  The  value  of  Kyfgnt  required  to  match  the  empiri¬ 
cal  lift  loss  in  ground  effect  is  shown  in  Figure  16  as  a  function 
of  plate  diameter,  D/De.  A  value  of  Kwent  =  1.50  was  selected  for 
use  in  all  future  cases  in  this  study.  Figure  17  presents  a 
summary  of  the  procedure  used  to  calculate  the  total  suckdown 
force  on  V/STOL  aircraft  hovering  in  ground  effect. 


2.  FOUNTAIN  FORCE  CALCULATIONS  -  Using  the  MCAIR  V/STOL  force 
prediction  methodology,  the  total  multiple  jet  fountain  upwash 
momentum  flux  incident  upon  the  aircraft,  rtf  j ,  is  calculated  as 
the  sum  of  an  inner  region  contribution  rtfj^n  and  an  outer  region 

contribution,  ^flout*  The  fountain  upwash  momentum  flux 
contributions  emanating  from  the  inner  region  and  outer  region  are 
determined  graphically  (Figure  18)  using  the  data  determined  using 
the  AGSF  methodology.  However,  only  a  fraction  of  the  total 
impinging  fountain  momentum  flux  is  converted  to  a  force  on  the 
aircraft.  This  fraction  is  accounted  for  in  the  methodology  by 
employing  the  fountain  impingement  momentum  flux  coefficient,  £t. 
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Figure  14.  Paneled  Representation  ot  Axlsymmetric  Jet  Issuing  from  a  Circular  Plate 

in  Ground  Effect 
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Figure  15.  Wall  Jet  Entrainment  Velocity  Correction  Factor 
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Figure  16.  Wall  Jet  Entrainment  Velocity  Correction  Factor  Required 
to  Match  Empirical  Lift  Loss 

Circular  Plate  IGE 


MCAIR  SUCKDOWN  METHODOLOGY: 

1.  PANEL  AIRCRAFT  AND  LIFT  JETS 

Z  DETERMINE  (AL/T)CALC  OUT  OF  GROUND  EFFECT  (OGE)  USING 
UNCORRECTED  VpENT  OF  KLEIS  AND  FOSS 

3.  CORRECT  VpENT  TO  MATCH  EMPIRICAL  LIFT  LOSS,  (AL/T)EXp,  AS: 

„  /|AL/T)Exp 

Vfent'  Vs^wVfent 

4.  PANEL  WALL  JETS 

5.  DETERMINE  VyyENT  USING  AGSF  PROGRAM  AND  CORRECT  AS: 

VWENT*  1-SVwENT 

&  CALCULATE  TOTAL  SUCKDOWN  FORCE  IGE,  FS,  USING  PANEL 
METHOD  WITH  VpENT  AND  V^ENT 

QfMS4t«n« 

Figure  17.  Summary  of  Procedure  for  Calculation  of  Sucfcdown  Force  in  Ground  Effect 
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FOUNTAIN  IMPINGEMENT  LINES 


a.  Inner  Region  Fountain  Impingement  Momentum  Flux  -  To  cal¬ 
culate  *f“  ",  all  of  the  wall  jet  momentum  flux  within  the  inner 

region  is  assumed  to  be  turned  upward  in  a  central  fountain.  Due 
to  the  lack  of  sufficient  experimental  data  to  define  multiple  jet 
central  fountain  characteristics,  the  central  fountain  is  assumed 
to  be  circular  in  cross  section  and  to  e.inibit  the  same  spreading 
characteristics  as  described  below  for  the  two- jet  fountain.  The 
total  momentum  flux  contained  within  the  central  fountain  upwash 
is  determined  as  (see  Figure  19 ); 


«f 


in 


n  A04 

il  v  "in 


m 


where  Mje  is  the  jet  exit  momentum  flux,  X  ^  is  the  fountain 
formation  momentum  flux  recovery  factor  determined  for  two-jet 
fountains  in  Reference  5  and  n  is  the  number  of  lift  jets.  (The 
factor  XjCj  was  found  to  be  equal  to  approximately  .55  for  two  jets 
impinging  vertically  with  equal  thrust  and  this  value  increased  as 
the  thrust  differential  between  the  two  impinging  jets  increased; 
see  Figure  2.)  The  inner  region  fountain  impingement  momentum 
flux  is  then  calculated  based  on  the  ratio  of  the  cross-sectional 
area  of  the  central  fountain  at  impingement,  Afc,  and  the  planform 
area  intercepted  by  the  central  fountain,  Aff  as  (see  Figure  20): 
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Figure  19.  Inner  Region  Fountain  Momentum  Flux 
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gun  20.  Fountain  Impingement  Forces 

Uniform  Momentum  Flux  Distribution 


NADC-81106-60 


b.  Outer  Region  Fountain  Impingement  Momentum  Flux  -  To  cal¬ 
culate  ftflout'  the  fountain  upwash  momentum  vectors  are  projected 
upward  from  the  ground  plane  stagnation  line  to  intersect  the  air¬ 
craft  planform  yielding  the  fountain  impingement  centerline 
(Figure  18).  Application  of  an  enpirical  two-jet  fountain  spread¬ 
ing  rate  yields  the  effective  fountain  upwash  boundaries  and, 
consequently,  the  fountain  "footprint"  or  that  area  of  the  air¬ 
craft  planform  intercepted  by  the  fountain  upwash.  The  two- jet 
fountain  spreading  rate  was  determined  based  on  the  data  obtained 
in  Reference  5  (see  Figure  21)  and  is  given  as: 


where 


D 


.50  +  .35 


1 

D 


(9) 

(10) 


Previously,  the  fountain  impingement  momentum  flux  in  the  outer 
region  was  also  calculated  using  the  simple  area  ratio  between  the 
fountain  “footprint"  area  Aff,  and  the  fountain  cross-sectional 
area  at  impingement,  Afc,  (Figure  20)  with: 

B,_  Aff  B,  .... 

fIout  x  fnout 

where  Mfnout  represents  the  vertical  momentum  flux  contained  with¬ 
in  the  fountain  upwash  in  the  outer  region.  This  definition  of 
Bf^out  inherently  assumes  a  uniform  distribution  of  momentum  flux 
across  the  fountain.  The  apparent  overlap  of  the  cross-hatched 
areas  of  the  inner  region  and  outer  region  computational  models 
shown  in  Figure  20  accounts  for  the  mixing  between  the  two-jet 
upwash  flow  and  the  central  fountain  flow. 

Recently,  however,  a  more  accurate,  non-uniform  fountain 
momentum  flux  distribution  curve  was  developed  under  IRAD  based  on 
a  correlation  of  existing  2-jet  fountain  upwash  dynamic  pressure 
profiles  obtained  in  Reference  5  and  in  an  unpublished  IRAD  experi¬ 
mental  program.  Figure  22  presents  the  dynamic  pressure  profiles 
obtained  foij  a  nozzle  exit  spacing  of  s/D  *  .12.8  and  nozzle  exit 
height  of  h/D  *  5.0  and  for  s/D  *  7.49  and  h/D  *  8.0.  The  figure 
illustrates  that  the  dynamic  pressure  profiles  are  similar  for  var¬ 
iations  in  fountain  height  ^  and,  in  fact,  are  essentially  similar 
for  variations  in  s/D  and  h/D  as  well.  Consequently,  based  upon 
these  correlations,  a  single,  empirical  2- jet  fountain  dynamic 
pressure  profile  (Figure  23)  was  derived.  This  empirical  fountain 
dynamic  pressure  profile  was  then  integrated  with  respect  to  the 
fountain  width,  n*  to  yield  the  2- jet  fountain  upwash  momentum 
flux  distribution  shown  in  Figure  24.  The  decrease  in  slope  of 
the  distribution  curve  shown  in  Figure  24  with  n  indicates  the 
decreased  impact  of  the  fountain  upwash  along  the  fountain  outer 
edges,  a  trend  not  previously  accounted  for  with  the  uniform 
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fountain  momentum  flux  distribution.  This  2-jet  fountain  upwash 

momentum  flux  distribution  curve  was  used  to  graphically  determine 

ftfl  *8: 

“out 


fnOUtl 


ftf  (n) 

Af.  .  (nf}  nout 

tot 


Where  &fnout  is  the  fountain  normal  momentum  flux  incident  upon 

the  aircraft  in  the  outer  region  and  - -  was  obtained  from 

Figure  24.  ”^tot 

The  total  fountain  impingement  momentum  flux  was  then  calculated 


Afl  -  Aflin  +  L*flc 


v 


BASED  ON  EXPERIMENTAL  DATA  EFFECTIVE  FOUNTAIN  WIDTH 

OBTAINED  UNOER  CONTRACT  TO  ONR:  DEFINED  AS: 

V2  *  0.50  +  0.35  £  _  V2 

D  0  D  '  2  D 
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Figure  21.  2-Jet  Fountain  Spreading  Characteristics 
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Fig un  22.  2-Jat  Fountain  Dynamic  Praaaura  Profiles 
Equal  Thrust 


Figure  23.  Empirical  2-Jat  Fountain  Dynamic  Preasura  Proflla 


Figure  24.  2-Jat  Fountain  Upwash  Momentum  Flux  Distribution 
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SECTION  III 
CALCULATION  OF  X f j 


Using  the  approach  described  in  Section  II,  the  fountain 
impingement  momentum  flu*  transfer  coefficient,  XfI,  was  calcu¬ 
lated  for  the  YAV-8B,  the  MCAIR  Model  260  and  the  Grumman  Mo^el 
698  V/STOL  aircraft  at  three  aircraft  heights  above  ground, 

The  height,  h,  is  measured  normal  to  the  ground  plane  surface  £o 
the  point  on  the  fuselage  undersurface  directly  beneath  the  air¬ 
craft  center  of  gravity.  The  diameter,  De,  is  the  equivalent 
single  nozzle  exit  diameter  defined  in  equation  (6). 

The  primary  quantities  calculated  using  the  MCAIR  V/STOL 
force  prediction  methodology,  as  well  as  the  experimentally 
induced  lift  determined  for  each  of  the  aircraft  configurations, 
are  presented  in  Table  1.  The  suckdown  forces,  fountain 
impingement  momentum  flux  magnitudes  and  XfIs  calculated  for  the 
three  aircraft  are  discussed  below. 

1.  YAV-8B  -  The  fountain  impingement  momentum  flux  transfer  coef¬ 
ficient  was  previously  determined  for  the  YAV-8B  in  hover  IGE  at 
an  aircraft  height  of  h/D~  ■  2.76  in  a  related  MCAIR  Independent 
Research  and  Development  (IRAD)  study.  This  study  was  extended  in 
the  current  program  to  include  calculations  of  X  f  j  for  the  YAV-8B 
at  h/De  »  1.72,  3.80  and  4.84.  Figure  25  presents  the  paneled 

representation  of  the  YAV-8B  in  the  clean  configuration  without 
wing  pylons,  gun  pods  or  fuselage  strakes.  The  horizontal  and 
vertical  tail  surfaces  were  assumed  to  contribute  little  to  the 
overall  suckdown  force,  and  these  surfaces  were  omitted  from  the 
paneled  model  to  reduce  the  paneling  effort  and  computational 
costs.  This  groundrule  was  adopted  for  the  paneling  of  all  three 
aircraft.  For  the  in-ground-effect  analyses,  the  wall  jet  sur¬ 
faces  were  represented  as  flat,  zero- thickness  plates  extending 
approximately  60  nozzle  exit  diameters  downstream  from  the  jet 
impingement  points. 

Unfortunately,  no  experimental  data  are  available  for  the 
YAV-8B  aircraft  in  a  clean  configuration  (without  gun  pods  or 
LID's).  The  fuselage  contours,  nozzle  sizes  and  nozzle  locations, 
however,  are  virtually  identical  on  both  the  AV-8A  and  the  YAV-8B 
aircraft.  Consequently,  the  induced  lift  on  the  two  aircraft  was 
assumed  to  be  the  same,  and  the  AV-8A  data  were  utilized  for  the 
calculations  of  X^j  for  the  YAV-8B  in  the  clean  configuration. 
Figure  26  presents  the  experimental  data  obtained  for  the  AV-8A 
(from  Reference  13)  in  the  clean  configuration  and  for  the  YAV-8B 
(from  Reference  14)  configured  with  two  16  inch  fuselage  strakes. 
Also  shown  in  Figure  26  is  the  variation  with  height  of  the  suck¬ 
down  force,  Fg/T ,  and  the  incident  fountain  momentum  flux,  rt^j/T, 
calculated  using  the  MCAIR  V/STOL  methodology.  Figure  26  indi¬ 
cates  that  approximately  19%  of  the  total  jet  exit  momentum  flux 
(thrust)  is  recovered  in  the  fountain  upwash  at  h/De  *  1.72  and 
that  this  value  decreases  to  13%  at  h/De  3  4.84. 
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TABLE  1.  CALCULATED  FORCE  QUANTITIES 


Out  of  Ground  Effect 

Aircraft 

/AL\ 

' T  /exp 

fa 

T 

V'fant 

Vfent 

YAV-8B 

-0.0200 

-0.0044 

2.132 

Model  260 

-00080* 

-0.0033 

1.633 

Model  698 

-0.0101* 

-0.0003 

5.491 

'From  Byrne’s  correlation 


Mfl 

T 


In  Ground  Effect 


Aircraft 

Configuration 

H 

nil 

n 

PH 

Mf,  . 

,  out 

Vl,l 

g 

u 

mm 

E9 

T 

T 

1.72 

-0.006 

-0.061 

0.129 

Esa 

0.291 

YAV-8B 

2.76 

-0.047 

0.125 

0.050 

0.234 

without  LIDs 

3.80 

-0.018 

-0.039 

0.107 

0.041 

0.148 

0.142 

4.84 

-0.027 

-0.034 

0.098 

0.032 

0.1 3C 

0.054 

1.72 

0.066 

-0.061 

0.129 

0.060 

0.189 

0.672 

YAV-8B 

2.76 

0.060 

-0.047 

0.125 

0.050 

0.175 

0.611 

with  LIDs 

3.80 

0.039 

-0.039 

0.107 

0.041 

0.148 

0.527 

4.84 

0.012 

-0.034 

0.098 

0.032 

0.130 

0.354 

MrvlAl  2fi n 

0.92 

0 

-0.048 

0.064 

0.031 

0.095 

0.505 

2.42 

0.010 

-0.010 

0.040 

|'X'  <1 

0.063 

0.317 

3.92 

-0.009 

-0.010 

0.028 

0.044 

0.023 

Model  260 

0.92 

0.033 

-0.048 

0.064 

0.031 

0.095 

0.853 

with  LID$ 

2.42 

0.014 

-0.010 

0.040 

0.023 

0.063 

0.381 

3.92 

-0.007 

-  0.010 

0.028 

0.016 

0.044 

0.068 

0.48 

-0.027 

-0.037 

— 

0.101 

0.101 

0.099 

Model  698 

0.83 

-0.015 

-0.028 

— 

0.089 

0.089 

0.146 

with  Round 

1.89 

-0.004 

-0.017 

— 

0.064 

0.064 

0.203 

Chines 

2.50 

-0.003 

-0.015 

— 

0.054 

0.054 

0.222 

3.81 

-0.005 

-0.012 

— 

0.040 

0.040 

0.175 

0.48 

-0.004 

-  0.037 

— 

0.101 

0.101 

0.327 

Model  698 

0.63 

0.071 

-  0.028 

— 

0.089 

0.089 

1.112 

with  Square 

1.89 

0.025 

-0.017 

— 

0.064 

0.064 

0.656 

Chines 

2.50 

0.009 

-0.015 

— 

0.054 

0.054 

0.444 

3.81 

-0.005 

M'iK 

— 

0.040 

0.040 

0.175 

anurm-n 
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The  fountain  impingement  momentum  flux  coefficient  was  cal¬ 
culated  for  the  YAV-8B  by  rearranging  equation  (1)  to  yield: 


xfl 


AL/T  -  Fa/T 

Afl/T 


(14) 


Figure  27  presents  the  variation  of  Xfj  with  h/De  for  the  YAV-8B 
in  both  the  clean  configuration  and  with  the  16  inch  fuselage 
strakes.  The  presence  of  the  strakes  was  assumed  to  have  negli¬ 
gible  effect  on  the  calculated  values  of  Fg/T  and  fifi/T. 
Consequently,  the  curves  presented  in  Figure  26  were  used  to 
calculate  Xfj  for  both  configurations.  Figure  27  indicates  that 
only  29%  of  the  momentum  flux  within  the  fountain  upwash  is  con¬ 
verted  to  a  fountain  impingement  force  on  the  clean  configured 
aircraft  at  h/De  =  1.72.  At  h/De  **  4.84,  the  fountain  effects  are 
shown  to  be  practically  nonexistent  on  the  clean  aircraft  with 
a  0.05.  The  effect  of  adding  the  fuselage  strakes  is  seen  to 
increase  the  effectiveness  of  the  fountain  impingement  process  by 
capturing  additionally  approximately  38%  of  the  fountain  upwash 
momentum  flux. 


Figure  27.  Fountain  Impingement  Momentum  Flux 
Transfer  Coefficient 
YAV-8B 
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2.  MCAIR  MODEL  260  -  The  fountain  impingement  momentum  flux 

transfer  coefficient  was  determined  for  the  Model  260  aircraft  at 
aircraft  heights,  h/De,  of  .92,  2.42  and  3.92.  The  calculations 
were  made  for  the  aircraft  in  both  a  clean  configuration  and  in  a 
configuration  with  lift  improvement  devices  (LIDs)  consisting  of 
two  vertical  longitudinal  fuselage  strakes  and  a  forward  lateral 
vertical  fence.  Figure  28  presents  the  paneled  representation  of 
the  Model  260  aircraft  (without  tail  surfaces)  in  the  clean,  hover 
configuration  that  was  utilized  for  the  suckdown  calculations. 
For  the  IGE  analyses,  the  wall  jets  were  again  paneled  as  zero 
thickness  flat  plates  and  extended  approximately  30  nozzle  exit 
diameters  downstream  of  the  jet  impingement  points. 

The  experimentally  determined  induced  lift  (AL/T)  obtained 
using  a  4.1%  scale  model  of  the  Model  260  aircraft  (Reference  15), 
with  and  without  LIDs,  is  presented  in  Figure  29.  The  data  are 
shown  for  the  case  of  all  three  jets  impinging  vertically  with 
each  jet  operating  at  a  nozzle  pressure  ratio  of  1.50.  Also  shown 
in  the  figure  are  the  calculated  values  of  Fg/T  and  rtfj/T.  Figure 
29  indicates  that  the  Model  260  experiences  approximately  one 
third  of  the  suckdown  force  as  the  YAV-8B  and,  likewise,  recovers 
only  about  one  third  of  the  jet  exit  momentum  flux  in  the  fountain 
upwash . 

The  fountain  impingement  momentum  flux  transfer  coefficients 
calculated  for  the  Model  260,  both  with  and  without  LIDs,  are  pre¬ 
sented  in  Figure  30.  Once  again,  is  found  to  decrease  with 

increasing  h/D  .  For  the  configuration  without  LIDs,  approxi¬ 
mately  50%  of  the  impinging  fountain  momentum  flux  is  converted  to 
a  fountain  impingement  force  at  h/De  *  .92.  However,  at  h/De  = 
3.92  only  2%  of  the  impinging  fountain  upwash  momentum  flux  is 
effective  as  a  fountain  impingement  force.  The  LIDs  are  seen  to 
increase  the  efficiency  of  the  fountain  impingement  process  up  to 
an  aircraft  height  of  approximately  2.5Dg.  Above  this  height  the 
increase  in  is  less  than  .05  over  that  found  for  the  clean 

configuration. 

3.  GRUMMAN  MODEL  698  -  Calculations  of  were  made  for  the 

Grumman  Model  698  aircraft  at  aircraft  heights  of  h/De  =  .48,  .83 

1.89,  2.50  and  3.61.  The  effect  of  fuselage  corner  radius  was 
a  -so  investigated  by  examining  the  Model  698  with  both  round  and 
square  fuselage  corners  (or  chines).  Figure  31  presents  the 
paneled  representation  of  the  Model  698  in  the  cisan  hover 
configuration  (without  vertical  tail)  for  a  nozzle  exit  centerline 
spacing  (s/D)  of  2.8.  The  aircraft  nozzles  were  paneled  to 
simulate  the  experimental  test  installation. 

The  experimentally  induced  lift  obtained  from  Reference  16 
for  the  Model  690  with  both  round  and  square  chines  is  presented 
in  Figure  32.  The  data  was  obtained  for  the  4.2%  scale  model  with 
vertically  impinging  jets  operating  at  a  nozzle  pressure  ratio  of 
1.06.  Also  shown  in  Figure  32  are  the  suckdown  force,  Fg/T,  and 
the  fountain  momentum  flux  at  impingement,  flfj/T,  calculated  for 
the  Model  698. 
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Figure  28.  Paneled  Representation  of  MCAIR  Model  260 
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Figure  30.  Fountain  Impingement  Momentum  Flux 
Transfer  Coefficient 

Model  260 
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Figure  32.  Grumman  Model  698  Induced  Lift  in  Ground  Effect 

The  fountain  impingement  momentum  flux  transfer  coefficients 
calculated  for  the  Grumman  Model  698  aircraft,  with  both  round  and 
square  fuselage  chines  are  shown  in  Figure  33.  The  effect  of 
fuselage  corner  radius  was  assumed  to  have  negligible  effect  on 
the  suckdown.  Consequently,  the  values  of  Fp/T  and  Mfj/T  shown  in 
Figure  32  were  utilized  for  both  configurations  to  calculate 
For  the  aircraft  with  the  round  chines,  Xfj  displays  an  uncharac¬ 
teristic  decrease  with  decreased  h/De  below  h/D.a  ■  2.5.  For  the 
aircraft  with  the  square  chines,  A  f r  displays  the  expected  trend 
of  increasing  with  decreasing  h/De  down  to  an  aircraft  height  of 
h/De  *  -83.  Below  this  height,  however,  A  *  j  displays  an  unrealis¬ 
tic  decrease  in  magnitude  and  is  therefore  questionable.  The 
advantage  of  square  chines  over  round  chines  is  evident  from  both 
Figure  32  and  Figure  33.  A  possible  explanation  for  this  effect 
is  given  by  Kalemaris  in  Reference  16. 


31 


NADC-81106-60 


Figure  33.  Fountain  Impingement  Momentum  Flux  Transfer  Coefficient 

Grumman  Model  698 
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SECTION  IV 
SUMMARY  RESULTS 


The  fountain  impingement  momentum  flux  transfer  coefficients 
calculated  for  the  YAV-8B,  the  MCAIR  Model  260  and  the  Grumman 
Model  698  aircraft  are  summarized  in  Figure  34.  Since  the  present 
methodology  assumes  distinct  regions  of  free  jet  flow,  jet  impinge¬ 
ment,  wall  jets,  fountain  formation  and  fountain  upwash;  predic¬ 
tion  of  If i  for  h/De  <  1  have  been  omitted  from  Figure  34  since 
distinct  flow  regimes  do  not  exist  at  these  low  heights.  The 
sensitivity  of  Afj  due  to  fuselage  undersurface  shape  is  apparent 
by  comparing  the  curves  in  Figure  34  for  the  YAV-8B  with  and 
without  LIDs.  This  sensitivity  is  also  displayed  by  the  Model  260 
and  the  Model  698  at  the  lower  aircraft  heights  where  significant 
fountain  capture  exists . 

The  fuselage  cross-sections  associated  with  the  Model  260  and 
the  Model  698  aircraft  are  in  general  quite  similar,  displaying  a 
flat  bottom  with  rounded  corners.  The  YAV-8B  fuselage,  on  the 
other  hand  is  rounded  across  the  entire  width.  However,  the  foun¬ 
tain  impingement  momentum  flux  coefficients  calculated  for  the 
Model  260  (without  LID's)  are  quite  similar  to  those  determined 
for  the  YAV-8B  (without  LID's),  while  the  Xfjs  calculated  for  the 
Model  698  (with  round  chines)  were  much  lower.  This  suggests  that 
the  Xfx  associated  with  multiple  jet  V/STOL  aircraft  producing  a 
central  fountain  may  be  fundamentally  different  from  those  asso¬ 
ciated  with  2- jet  V/STOL  aircraft  where  no  central  fountain  is 
formed . 

Figure  34  also  illustrates  the  dependence  of  Xfj  on  h/De, 
with  Xfx  decreasing  with  increased  aircraft  height.  \f j  was 
originally  expected  to  be  relatively  insensitive  to  variations  in 
h/De  providing  the  fountain  upwash  was  adequately  modeled.  The 
h/De  dependency  illustrated  in  Figure  34  may  thus  be  a  result  of 
inadequate  fountain  modeling.  This  inaccuracy  may  result  from  one 
of  two  identified  problem  areas:  (1)  insufficient  empirical  data 
defining  central  fountain  core  upwash  characteristics,  and/or  (2) 
lack  of  an  adequate  model  defining  conditions  for,  and  magnitude 
of,  fountain  upwash  interaction  with  the  downward  impinging  jet 
flows.  (The  characteristics  of  four- jet  central  fountains  were 
investigated  in  a  MCAIR  Independent  Research  and  Development 
experimental  program  in  1982. )  The  utility  of  the  MCAIR  V/STOL 
force  prediction  methodology,  however,  is  not  severely  impaired  by 
the  absence  of  these  fountain  models  since  inadequacies  associated 
with  the  fountain  upwash  models  are  inherently  accounted  for  in 
Xfj. 
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Figure  34.  Fountain  Impingement  Momentum  Flux  Transfer  Coefficient  Summary 
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The  agreement  between  the  indirect  and  direct  methods  of 
measuring  mass  entrainment  discussed  in  the  Appendix  and  shown  in 
Figure  A13  indicates  that  both  approaches  are  accurate.  Although 
the  entrainment  rates  determined  in  this  study  are  higher  than 
those  determined  previously  (Figure  Al),  the  difference  is  still 
too  small  to  account  for  the  magnitude  of  the  suckdown  forces 
reported  in  Section  II. 1. a*. 


*The  suckdown  force  on  the  Model  698  aircraft  out-of-ground  effect 
was  recalculated  using  the  entrainment  rates  determined  for  nozzle 

A  at  NPR  *  1.10  (m^  *  .334,  m2  *  .422)  as  ^  *  -.0008.  This  may 

be  compared  with  the  semi-empirical  value  of  -t*  ■  -.0101 

determined  using  Byrne's  formula  (Figure  10). 
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SECTION  V 
CONCLUSIONS 


Based  on  the  calculations  of  the  suckdown  force  determined 
for  the  MCAIR  YAV-8B,  the  MCAIR  Model  260  and  the  Grumman  Model 
698  in  hover  out  of  ground  effect,  it  was  concluded  that  existing 
empirical  free  jet  entrainment  rates  are  too  weak.  The  results  of 
the  experimental  investigation  of  free  jet  mass  entrainment  vali¬ 
date  the  indirect  method  and  indicate  that  the  indirect  method  may 
be  applied  to  situations  in  which  use  of  the  direct  method  may  be 
impractical.  Three  such  situations  are  depicted  in  Figure  35. 
The  first  situation  represents  the  case  of  a  high  velocity  jet  in 
which  the  flow  dynamics  would  be  too  severe  for  penetration  with  a 
standard  hot  film  probe.  The  second  situation  is  that  of  two  or 
more  merging  jets.  In  this  case,  distinct  free  jet  velocity 
profiles  may  not  be  attainable  or  the  required  survey  volume  may 
be  too  large  for  use  of  the  direct  method.  The  final  situation 
represents  the  case  of  a  very  hot  jet  in  which  immersion  by 
standard  probes  may  not  be  practical . 


The  Xfjs  determined  for  the  three  aircraft  using  the  MCAIR 
hover  methodology  with  adjusted  free  jet  and  wall  jet  mass  entrain¬ 
ment  rates  indicated  a  general  trend,  with  \fj  decreasing  with 
increased  aircraft  height  above  ground.  Use  of  additional  lift 
jets  decreases  the  rate  of  decay  of  Xfj  with  aircraft  height. 


Jets 


(Mj  >  1.0) 


Jets 

Tj  >  >  Tamb 


Figure  35.  Potential  Appficatlont  of  the  Indirect  Method  of  Determining  Free-Jet 

Maes  Entrainment 
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APPENDIX  A 

During  the  study  described  in  the  preceding  sections,  it  was 
determined  that  existing  empirical  axisymmetric  free  jet  entrain¬ 
ment  rates  (Figure  Al)  are  too  low  for  the  accurate  prediction  of 
suckdown  forces  on  V/STOL  aircraft  in  hover.  The  conventional 
experimental  approach  for  determining  the  entrained  mass  distri¬ 
bution  is  based  upon  detailed  velocity  measurements  within  the 
jet.  If  the  axial  velocity  component  is  measured  at  enough  points 
on  a  cross-section  to  establish  a  velocity  profile,  then  the  local 
mass  flux  can  be  determined  by  integration.  The  accuracy  of  this 
direct  approach,  however,  is  highly  sensitive  to  the  accuracy  of 
the  probe  sensor  and  survey  position  calibration  and  to  the 
numerical  integration  method  selected.  As  an  alternative  to  this 
approach,  MCAIR  developed  a  new  technique  to  indirectly  measure 
the  entrainment  rates  associated  with  axisymmetric  free  jets. 
This  indirect  method  was  developed  as  a  means  of  validating  the 
direct  approach  and  to  verify  the  potential  flow  model  of  jet 
entrainment  induced  flowfields.  The  indirect  method  utilizes  a 
hot  film  anemometer  (HFA)  to  measure  the  induced  velocity  field 
about  subsonic,  axisymmetric  free  jets.  Inverse  potential  flow 
theory  is  employed  to  convert  the  measured  induced  velocities  to 
mass  entrainment  rates.  This  technique  was  applied  in  an  experi¬ 
mental  program  by  MCAIR  in  order  to  quantify  the  dependence  of 
free  jet  mass  entrainment  on  nozzle  pressure  ratio  and  turbulence 
intensity  for  isolated  subsonic  axisymmetric  jets.  The  technical 
approach,  experimental  test  program  and  test  results  are  described 
in  the  following  sections. 


Entrained 
Mass  Flow 


Axial  Distance  from  Nozzle  P'ane  (z  /D) 


i 


9nxnM 


Figure  Al.  Entrained  Maes  Flow  Distribution  lor  Circular  Jets 
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A.l  TECHNICAL  APPROACH  -  Consider  a  set  of  n  points  outside  of 
the  mixing  region  of  a  jet  (Figure  A2).  Suppose  that  at  each 
point  (zi»ri),  the  radial  velocity  component  Vr  has  been  measured 
experimentally.  The  data  reduction  procedure  described  below  will 
determine  the  jet  entrained  mass  distribution  that  most  nearly 
corresponds  to  the  measured  velocity  field,  in  the  sense  of  a 
least  square  error. 


Figure  A2.  Representation  ot  an  Axisymmetrlc  Subsonic  Jet  for  Indirect 
Mass  Entrainment  Measurements 


The  mechanism  for  converting  measured  velocity  to  mass 
entrainment  is  inverse  potential  flow  theory.  It  is  assumed  that 
the  jet  can  be  represented  by  a  potential  core  region  of  length  zc 
and  a  semi- infinite  fully  developed  region  (Figure  A2).  Prior 
investigations  have  shown  zc  to  be  approximately  equal  to  six 
nozzle  exit  diameters  (zc  =  6D).  The  rate  of  mass  entrainment  per 
unit  length  (dm/dz)  is  assumed  to  be  an  unknown  constant  nq  in  the 
potential  core  region  and  a  second  unknown  constant  m2  in  the 
fully  developed  region.  If  the  values  and  m2  were  known  a 
priori,  the  calculation  of  induced  radial  velocity  component  Vr  at 
any  point  (z,r)  would  be  straightforward.  The  jet  would  be 
represented  by  a  distribution  of  sinks  along  the  centerline,  where 
the  sink  strength  would  be  mL  in  the  potential  core  region 
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(0_<z<zc)  and  m2  in  th©  fully  developed  region.  Image  sinks  for 
z<0  would  simulate  the  presence  of  the  flat  exit  plane  at  z  =  0. 
The  formula  for  Vr  is: 


Vr  (z,r) 


m2 

2irrp 


mi  -  m2 

- 

(z  +  ZC) 

(z  -  zc) 

4  TOP 

/  r*  +  (z  +  z  ) ^ 

/  r^  +  (z  -  zc) ^ 

(Al) 


where  the  density  p  is  constant  in  low  speed  flow. 


Here  it  is  desired  to  solve  for  the  unknowns  mi  and  m2  such 
that  Vr  (z,r)  from  equation  (Al)  will  approximately  match  the 
experimentally  measured  value  at  each  of  the  probe  positions 

(zi^i).  Formally,  the  objective  is  to  minimize  the  quadratic 
error  function. 


E 


n 

-  Z 
i=i 


Vr(zi,ri)  - 


2 


(A2) 


The  minimization  of  E  is  accomplished  by  substituting  (Al) 
into  (A2)  and  differentiating  (A2)  with  respect  to  the  two 
unknowns.  This  results  in  a  system  of  two  equations  which  may  be 
solved  directly  for  mi  and  m2.  Following  the  determination  of  irq 
and  m2  the  entrained  mass  distribution  can  be  expressed  as 


i-miz 

m  (mi  -  m2 )  zc  +  nvjz 


CKz_<Zc 

zc<z 


(A3) 


A.  2  EXPERIMENTAL  TEST  PROGRAM  -  The  experimental  test  program  was 
conducted  in  the  Inlet  Simulator  test  cell  of  the  MCAIR  Propulsion 
Subsystem  Test  Facility  (PSTF) .  The  facility  utilizes  high 
pressure  (4137  kpa)  air  from  the  air  storage  system  of  the  MCAIR 
Polysonic  Wind  Tunnel.  The  experimental  program  was  designed  to 
provide  accurate  measurements  of  the  mass  entrainment  induced  by 
an  axisymmetric  free  jet  for  parametric  variations  in  nozzle 
pressure  ratio  (NPR)  and  free  jet  turbulence  intensity  levels. 
.The  experimental  test  hardware,  test  procedure  and  results  are 
described  below. 


A. 2.1  Test  Hardware  and  Instrumentation  -  Two  existing 
converging  nozzle  pipes  were  used  for  the  test  program.  These 
nozzles,  designated  nozzle  A  and  nozzle  B,  are  illustrated  in 
Figure  A3.  Nozzle  A,  which  was  used  for  the  investigation  of 
nozzle  pressure  ratio  effects,  has  a  nozzle  exit  diameter  (DA)  of 
2.87  cm.  Nozzle  B  was  used  for  the  investigation  of  turbulence 
effects  and  has  a  nozzle  exit  diameter  (Dg)  of  5.90  cm.  Each  of 
the  nozzles  is  instrumented  with  four  (4)  nozzle  exit  static 
pressure  taps,  one  (1)  total  pressure  probe  and  one  (1)  total 
temperature  probe.  In  addition,  nozzle  B  is  constructed  such  that 
a  series  of  turbulence  generating  plates  and  screens  may  be 
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Inserted  in  the  nozzle  pipe  at  a  distance  of  3.23  Dg  (19.05  cm) 
upstream  of  the  nozzle  exit  plane.  Substantial  effort  was  spent 
to  determine  the  most  suitable  combinations  of  plates  and  screens 
for  producing  turbulence  levels  different  from  that  of  the  base¬ 
line  nozzle  configuration  and,  at  the  same  time,  yielding  a  sym¬ 
metrical  jet.  The  final  nozzle  configurations  utilized  for  the 
turbulence  effects  investigation  are  designated  TO,  T7  and  T8. 
Nozzle  configuration  TO  is  the  baseline  configuration  and  consists 
of  nozzle  B  with  no  internal  plates  or  screens.  Nozzle  configura¬ 
tion  T7  consists  of  the  baseline  configuration  with  the  addition 
of  a  perforated  plate  upstream  of  the  nozzle  exit.  Nozzle  configu¬ 
ration  T8  is  represented  by  the  T7  nozzle  assembly  with  a  screen 
(T5 )  installed  directly  downstream  of  the  perforated  plate.  The 
T7  plate  and  the  T5  screen  utilized  for  configuration  T8  are 
illustrated  schematically  in  Figure  A4. 

A  127.0  cm  diameter  circular  nozzle  exit  plane  plate  was  fab¬ 
ricated  from  plexiglass  for  the  test.  Provisions  were  made  to 
enable  the  attachment  of  a  remotely  controlled  single  degree-of- 
freedom  probe  traverse  device  to  the  back  of  the  plate  such  that 
the  hot  film  anemometer  probe  could  be  traversed  in  a  direction 
normal  to  the  plate  surface  at  any  one  of  nine  radial  locations. 
The  nozzle  exit  plane  plate  was  adaptable  for  use  with  either 
nozzle  A  or  nozzle  P. 

Free  jet  mean  velocity  (U)  and  turbulence  intensity  (u'/U) 
levels  were  established  using  a  single  sensor  TSI  Model  1213-20 
hot  film  anemometer  probe.  The  jet  entrainment  induced  velocity 
field  surrounding  the  axisymmetric  jet  was  mapped  out  using  a 
specially  designed  TSI  Model  1330XX  temperature  compensated  hot 
film  probe  (Figure  A5).  In  order  to  maximize  the  sensitivity  of 
the  probe  to  the  radial  inflow,  the  hot  film  sensor  was  placed  on 
the  upstream  side  of  the  probe  only.  In  addition,  since  the  hot 
film  probe  is  primarily  sensitive  to  the  velocity  component  normal 
to  the  sensor  axis,  the  probe  was  oriented  such  that  the  sensor 
axis  was  parallel  to  the  free  jet  centerline.  The  hot  film  probes 
were  calibrated  using  a  TSI  Model  1125  calibrator  with  an  accuracy 
of  +1.0%  for  velocities  above  3.0  m/sec  and  with  an  accuracy  of 
+2.ff%  for  a  velocity  range  of  .15  m/sec  -  3.0  m/sec.  The  RMS 
value  of  the  fluctuating  portion  of  the  free  jet  velocities  (u') 
were  determined  using  a  Ballantine  Model  323  true  RMS  meter.  The 
free  jet  nozzle  exit  velocity  profiles  were  determined  using  a 
United  Sensors  Model  PFC-18-C-16-KL  pitot-static  pressure  probe. 
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Figure  A3.  Baseline  Nozzle  Configurations 


anMWt; 


m  a. 


NADC-811G  6-60 


Figure  A5.  TSI  Modal  1330XX  Temperature  Compensated  Hot  Film  Anemometer  Probe 


A. 2. 2  Test  Procedure  -  The  experimental  test  program  is  out¬ 
lined  in  Table  Al .  Nozzle  exit  velocity  profile  surveys  were  per¬ 
formed  to  investigate  the  uniformity  and  symmetry  of  the  nozzle 
exit  flow.  Figure  A6  presents  the  nozzle  exit  velocity  profiles 
obtained  for  nozzle  A  at  NPR  =  1.10,  1.35  and  1.80.  The  velocity 
U  is  normalized  by  the  jet  centerline  velocity  (Uql)  obtained  at 
y  =  0.  Figure  A6  illustrates  that  the  nozzle  exit  flow  was 
relatively  uniform  (within  2%  of  UCl)  for  all  three  nozzle 
pressure  ratios  except  near  the  edge  of  the  jet  where  viscous 
effects  occur.  The  effect  of  the  turbulence  generators  on  the 
velocity  profile  generated  by  nozzle  B  at  NPR  *  1.10  is 
illustrated  in  Figure  A7.  The  baseline  nozzle  (TO)  displays  a 
uniform,  symmetrical  "top  hat"  profile.  Addition  of  the  T7  plate 
accelerated  the  flow  at  the  edge  of  the  jet  to  approximately 
1.22  Uql*  Addition  of  the  T5  screen  downstream  of  the  T7  plate 
(configuration  T8)  smoothed  out  the  peaks  somewhat,  but  still 
produced  flow  at  the  edge  of  the  jet  that  was  approximately  8% 
higher  than  CJCL. 

u ' 

The  turbulence  intensity  (TT-)  of  the  jets  was  determined  for 
NPR  *  1.10  using  the  hot  film  anemometer  to  simultaneously  measure 
the  root  mean  square  of  _the  fluctuating  axial  velocity  ( u *  )  and 
the  mean  axial  velocity  (U).  (Extreme  flow  dynamics  prevented  use 
of  the  delicate  hot  film  anemometer  sensor  for  nozzle  A  at  the 
higher  NPRs . )  The  free  jet  mean  axial  velocity  profiles  and  turbu¬ 
lence  intensity  profiles  obtained  for  nozzle  B  in  the  TO,  T7  and 
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m8  configurations  at  NPR  *  1.10  and  2/DB  =  1.0  are  presented  in 
Figure  A8.  X-Y-Y  plotter  tracings  of  the  HFA  output  voltages  were 
used  to  interpolate  between  data  points  on  the  nean  velocity 
profile. 


TABLE  Al.  TEST  PROGRAM 


Exit  Velocity  Prof  lb  Surveys 

•  Nozzle  A 

-  NPR=1.10,  1.35,  1.80 

-  Z/Da  =  0.04 

—  Y/DA=-0.43,  -0.35,  -0.26,  -0.18,  -0.08,0,0.09,0.18,0.27,0.36,0.42 

•  Nozzle  8 

—  NPR  a  1.10 

—  Turbulence  Generator  TO,  T7,  T8 

-  Z/Db  =  0.02 

—  Y/Dga-0.47,  -0.43,  -0.35,  -0.26,  -0.17,  -0.08,0,0.09,0.17,0.26,0.35,0.44,0.47 

Turbulence  Intensity  Surveys 

•  Nozzle  A 

-  NPR=1.10 

-  Z/DA  *  1.0,  2.0,  4.0,  6.0,  9.0,  12.0 

•  Nozzle  B 

—  NPR  ss  1.10 

—  Turbulence  Generator  TO,  T7,  T8 

-  Z/Dg=  1.0,  2.0,  4.0,  8.0,  9.0 

Entrainment  Velocity  Measurements 

•  Nozzle  Pressure  Ratio  Effects 
—  Nozzle  A 

—  NPR  =1.10,  1.35,  1.80 

-  R/Da  =  4.0,  5.0,  6.0,  7.0,  8.0, 10.0, 12.0 

-  Z/Da  =  0.5,  1.0,  2.0,  3.0,  4.0,  6.0,  8.5,  11.0 

-  <*>  =  0‘,  180*,  270* 

•  Turbulence  Effects 
—  Nozzle  B 

-  NPR  =  1.10 

—  Turbulence  Generator  TO,  T7,  T8 

-  R/DB  =  2.913,  3.885,  4.856,  5.827,  7.284,  9.225 

-  Z/Dg  =  0.5,  1.0,  2.0,  3.0,  4.0,  5.5 

-  <*>  =  0*,  180*.  270* 

OnM7M.1t 
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The  free  jet  entrainment  induced  radial  velocities  were 
recorded  by  traversing  the  HPA  probe  in  the  axial  (z/D)  direction 
at  several  radial  (r/D)  stations.  Figure  A9  depicts  the  test  hard¬ 
ware  installation  used  to  measure  the  induced  velocities  for 
nozzle  A.  Data  points,  recorded  with  the  air  off  at  the  beginning 
and  end  of  each  axial  traverse  were  summed  and  averaged  to  provide 
a  correction  factor  to  account  for  any  extraneous  drafts  in  the 
test  cell.  This  average  air  off  velocity  (usually  on  the  order  of 
.0  3  n/sec)  was  subtracted  from  each  of  the  measured  induced 
velocities  recorded  during  the  traverse  to  yield  a  corrected 
entrainment  velocity.  Table  Al  defines  the  axial  and  radial 
stations  at  which  the  entrainment  velocities  were  recorded.  Due 
to  the  larger  diameter  of  nozzle  P  and  the  fixed  HFA  probe  length, 
surveys  of  the  nozzle  B  configurations  were  limited  to  Z/DB<  6. 

In  a  check  of  the  symmetry  of  the  entrainment  flow  about  the 
jet,  it  was  discovered  that  the  entrainment  velocities  measured  at 
an  azimuthal  station  of  $  *  0°  (see  Figure  A10)  were  generally 
lower  and  more  unsteady  than  those  measured  at  0  *  180°. 
Additional  measurements  made  at  0  =  270°  produced  results 
generally  falling  between  those  obtained  at  0  *  0°  and  0  =  180°. 
This  flow  asymmetry  was  found  to  persist  for  all  nozzle  configura¬ 
tions  and  all  nozzle  pressure  ratios.  In  addition,  rotating  the 
nozzle  180®  on  the  plenum  produced  no  change  in  the  azimuthal 
distribution  of  the  asymmetry.  It  was  thus  concluded  that  the 
flow  asymmetry  was  a  result  of  the  test  cell  structure  and 
internal  aerodynamics.  Consequently,  measurements  of  the 
entrainment  velocities  were  recorded  at  azimuthal  stations  of  0  = 
0°,  180®  and  270®  for  each  configuration  to  provide  a  larger,  more 
representative  data  base.  (Hardware  restrictions  prevented 
surveys  at  the  0  =  90®  location.) 

A. 2. 3  Test  Results  -  The  data  reduction  technique  described 
in  Section  A. 2.1  was  utilized  in  conjunction  with  a  potential  core 
length,  zc,  of  6D  to  calculate  m^  and  m2*  (A  sensitivity  study 
indicated  ni  and  relatively  independent  of  zc  for 
4  <_  z^D  ^8.)  The  calculated  values  of  m^  and  m2  determined  for 
nozzle  A  are  presented  in  Figure  All  as  a'*’  function  of  NPR.  Also 
shown  in  the  figure  is  the  value  of  mj_  determined  for  nozzle  B  in 
the  baseline  TO  configuration  at  NPR  =  1.10.  As  seen  in  the 
figure  both  entrainment  rates  m^  and  m2  are  higher  than  the  values 
determined  by  Kleis  and  Foss  which  are  currently  used  in  the  MCAIR 
v/STOL  methodology. 


49 


NADC-8 1104-60 


3P23078*  I 


Figure  A9.  Induced  Velocity  Measurement  Test  Hardware 

Nozzle  A 
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Vi* w  Looking  Upstntm 


Figure  A10.  Probe  Azimuthal  Locations 


Since  the  entrainment  measurements  for  the  turbulence  effects 
investigation  using  nozzle  B  were  limited  to  z/D  <_  6,  it  was  not 
possible  to  measure  m2  directly.  Thus,  based  on  the  data  of  Kleis 
and  Foss,  m2  was  defined  as 


m2 


(A4) 


The  effect  of  the  magnitude  of  m^  on  the  calculation  of  m1  was 
investigated  and  found  to  be  relatively  insignificant.  Figure  A12 
presents  the  effect  of  turbulence  intensity  on  Bi  for  NPR  =  1.10. 
The  turbulence  intensity  shown  is  the  value  at  the  jet  centerline 
at  z/D  »  1.0.  As  seen  _from  the  figure,  there  seems  to  be  little 
correlation  between  u'/U  and  mj_*  However,  it  is  noteworthy  that 
mi  is  again  higher  than  the  value  determined  by  Kleis  and  Foss  for 
all  four  cases. 

The  mass  entrainment  induced  by  jets  issuing  from  both  nozzle 
A  and  nozzle  B  at  NPR  ■  1.10  was  also  calculated  under  IRAD 
funding  using  the  direct  method  of  integrating  the  free  jet 
velocity  profiles.  The  velocity  profiles  were  obtained  during  the 
free  jet  surveys  using  the  hot  film  anemometer.  Figure  A13 
compares  the  mass  entrainments  determined  using  both  the  indirect 
and  direct  methods  for  nozzle  A  and  nozzle  B.  The  agreement 
between  the  two  methods  is  seen  to  be  good,  particularly  in  the 
potential  core  region.  Also  shown  in  the  figure  for  comparison  is 
the  mass  entrainment  determined  by  Kleis  and  Foss. 
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Figure  A12.  Effect  of  Turbulence  Intensity  on 
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Figure  A13.  Comparison  of  Indirect  and  Direct  Methods  for  Determining 

Free-Jet  Entrainment 
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